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^ at LHC energies. Relativistic terms in the production amplitude connected with the 

D ' relative motion of heavy quarks and the transformation law of the bound state wave 

^Hj' functions to the reference frame of moving charmonia are taken into account. For the 
gluon and quark propagators entering the amplitude we use a truncated expansion in 

^ , relative quark momenta up to the second order. Relativistic corrections to the quark 

I bound state wave functions are considered by means of the Breit-like potential. It 

' turns out that the examined effects decrease total nonrelativistic cross section more 

■ than two times and on 20 percents in the rapidity region of LHCb detector. 
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■ The high energy and luminosity of the LHC makes it possible to observe the pair charmo- 
^ , nium production processes and to measure the corresponding cross sections with sufficiently 



high accuracy. The result of the measurement of the pair J/^p meson production by the 
LHCb Collaboration was published in |l| and discussed many times on different workshops 
[2|. This process together with the charmonium and associated open charm production can 
be considered as a probe of the quarkonium production mechanism. According to non- 
relativistic QCD (NRQCD) |3| and collinear parton model 0^ the predictions of observed 
cross section in the leading order in can be obtained by the use of parton distribution 
functions and a set of local non-perturbative charmonium production color singlet and color 
octet matrix elements [7-ll|. In the proton-proton collisions, additional contributions from 



other mechanisms, such as the double parton scattering (DPS) or the intrinsic charm content 
of the proton to the total cross section are possible |12|-|l4|. The processes of quarkonium 
production in proton-proton interaction are generally described using the scale-dependent 
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parton density functions. They are calculated as functions of the Bjorken variable x at 
some factorization scale within the approach of the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) evolution equations. However, the double charmonium production in pp 
collisions at high energies can be sensitive to the details of the parton kinematics. There- 
fore, it is more appropriate to use the parton distributions unintegrated over the transverse 
momentum kt in the framework of the fcj-factorization 121 . There exists another source of 



theoretical uncertainty related with the pair charmonium production which gives essential 
modification of the cross sections. It is connected with the account of relative motion of 
heavy quarks forming the bound states. As was shown in [l5[ these relativistic corrections 
significantly change the cross section of the pair charmonium production in pp interaction 
obtained in non- relativistic approximation. The detailed investigation of this relativistic 

evi- 



mechanism for exclusive double charmonium production in e^e~ annihilation |16|-|18 
dently shows that it is impossible to obtain the reliable theoretical predictions for observed 
quantities without an account of relativistic corrections. Finally, the next to leading order 
QCD corrections to the production amplitudes also should be taken into account. 

The strategy of experimental investigations can be directed on the study of such physical 
reactions in pp collisions in which one of the described mechanisms of quarkonium production 
is dominant. Unfortunately, as we know at present time all enumerated mechanisms have 
important effect in the pair charmonium production and their contributions to the total 
cross section should be taken into account to obtain high accuracy theoretical result. 

In this work we continue the study of relativistic effects in the inclusive pair charmonium 
production by considering the process p+p — )■ rjc + rjc + X. Our calculation of the production 
cross section is performed on the basis of relativistic quark model used previously for the 
investigation of relativistic corrections to the other reaction p+p — )■ J / 'ip+J / in [HI . We 
work within the single-parton scattering (SPS) mechanism in which the basic contribution 
to the charmonium production is determined by the gluon-gluon fusion. The aim of the 
present study consists also in the analysis of some uncertainties regarding the choice of 
parton distribution functions (PDF). In spite of existing difficulties in the detecting of r]c 
meson pairs it is thought that in new run of the LHC this process will be studied more 
successfully. 

The differential cross section da for the inclusive double charmonium production in 
proton-proton interaction can be presented in the form of the convolution of partonic cross 
section da[gg — ?■ rjcrjc] with the parton distribution functions of the initial protons 043]: 

da[p + p r]c + 7]^ + X] = j dxidx2fg/pixi, /i)/g/p(x2, fx)da[gg -> r]cr]c], (1) 

where fg/p{x,fi) is the partonic distribution function for the gluon in the proton, xi,2 are 
the longitudinal momentum fractions of gluons. The cross section formula ([1]) contains the 
factorization of the long distance PDFs and the short distance gluon fusion cross section 
da[gg — )■ rjcrjc] with the factorization scale /z. Neglecting the proton mass and taking the 
cm. reference frame of initial protons with the beam along the 2;-axis we can present the 
gluon on mass-shell momenta ki^2 = a^i, 2-^(1, 0, 0, ±1). is the center-of-mass energy in 
proton-proton collision. 

In the quasipotential approach the double charmonium production amplitude for the 
basic parton subprocess g + g ^ rjc + rj^ can be expressed as a convolution of a pertur- 
bative production amplitude of two c-quark and c-antiquark pairs T{pi,P2', 5'i,Q'2) and the 
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FIG. 1: The typical diagrams (the set of 31 Feynman diagrams) of the leading order for g + g ^ 
r]c{J I'lp) + r]c{J/ip)- The others can be obtained by reversing the quark lines or interchanging the 
initial gluons. 



quasipotential wave functions of the final mesons [l7|, ll8| : 



M[gg VcVc]{ki,k2,P,Q) 



dp 

(2^3 J (27r) 



T(pi,P2;gi,g2), (2) 



where pi and p2 are four-momenta of c-quark and c-antiquark in the pair forming the first 
7]c particle, and q2 and qi are the appropriate four-momenta for quark and antiquark in the 
second meson rjc- They are defined in subsequent transformations in terms of total momenta 
P{Q) and relative momenta p{q) as follows: 



pi,2 = ^P±p, (pP)=0; qi,2 = lQ±q, (qQ) = 



(3) 



This expression describes the symmetrical escape of the c-quark and c-antiquark from the 
mass shell. In Eq. ([2]) we integrate over the relative three-momenta of quarks and antiquarks 
in the final state. The systematic account of all terms depending on the relative quark mo- 
menta p and g in ([1]) is important for increasing the accuracy of the calculation, p = Lp{0,p) 
and q = LQ(0,q) are the relative four- momenta obtained by the Lorentz transformation of 
four- vectors (0, p) and (0, q) to the reference frames moving with the four- momenta P and Q. 

The relativistic wave functions of the bound quarks, accounting for the transformation 
from the rest frame to the moving one with four momenta P and Q, are the following 
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where the hat symbol means a contraction of the four-vector with the Dirac gamma matrices; 
Vi = P/M, V2 = Q/M; e{p) = a/ m'^ + p^, m is c-quark mass, and M is rjc charmonium 
mass, M ^ 2m. 

The amphtude ([2]) is projected onto a color singlet state by replacing Vi{0)uk{0) with 
a projection operator of the form Vi{0)uk{0) = 75(1 + 7o)5ifc/2-\/6. The relativistic wave 
functions in Eq. (jlj) are equal to the product of the wave functions in the rest frame "^q" and 



the spin projection operators that are accurate at all orders in |p|/m [l7, 18|. Our derivation 
of relation (jl]) accounts for the transformation law of the bound state wave functions from 
the res t frame to the moving one with four momenta P and Q, which was obtained in 
The physical interpretation of the double charmonium production amplitude is the 
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following: we have a complicated transition of two heavy c-quark and c-antiquark produced 
in gg-fusion outside the mass shell and their subsequent evolution firstly on the mass shell 
(free Dirac bispinors) and then to the quark bound states. In the spin projectors (jl]) we have 
p^, 7^ M^/4 — so that we can consider these structures as a transition form factors 
for the heavy quarks from the free states to the bound states. 

In the leading order in the strong coupling constant a^, there are 39 Feynman diagrams 
contributing to the pair production of rjc mesons. They can be divided into two sets shown 
in Figs. [U and El respectively. Their total contribution to the production amphtude (|2]) can 
be presented in the following form: 



M[gg ^ VcVc]{h,k2,P,Q) = ^MttW^ ^ ^ [Tr 971 + 3A9}t] , 



(5) 



where we explicitly extracted the relativistic normalization factors a/2M of quasipotential 
wave functions. The construction and transformation of the production amplitudes is per- 
formed by means of the package FeynArts j2l| for the system Mathematica and Form j22| . 

The integrand term in ([5]) containing the trace of the amplitude Wl represents the con- 
tribution of 31 diagrams in Fig. [1] and equals up to the wave functions definitions {M the 
analogous expression in the case of pair J/ip production, which can be found in Ref. [l5|. The 
second integrand term in (jS]), coming from additional 8 diagrams in Fig. [2] for rjc production 
amplitude, has the form 
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(6) 



X 



'{k2—piY — m? ' {k2—p2Y — rn? 

where the Mandelstam variables for the gluonic subprocess gg — )■ rjcTjc are 

S={k^ + k2f = {P + Qf = XiX2S, 

t = {P- kif = {Q- k2f = M^ - xiVS{Po - |P| COS0) = 



(7) 
(8) 
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FIG. 2: The additional typical diagrams (the set of 8 Feynamn diagrams) of the leading order for 
g + g ^ r]c + r]c only. 



= - X1X2S + X2VS{Po + |P| COS0), 

U = {P- k2f = {Q- hf = - X2V^(Po + |P| COS( 

= - X1X2S + 2;iv^(Po - |P| COS0), 



(9) 



is the angle between P and the 2;-axis. The transverse momentum Pt of rjc and its energy 



Pn can be written as 
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In order to calculate relativistic corrections contained in the production amplitude ([6]) we 
expand the inverse denominators of gluon and quark propagators as series in relative quark 
momenta p and q: 
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There are 16 different propagators in the amplitude (j5]), which have to be expanded in the 
manner of Eqs. (ITT!) . Then, preserving in the expanded amplitude terms up to the second 
order in relative quark momenta p and q, we can perform angular integration using the 
following relations for iS-wave charmonium: 



^o(p) dp 



I m 2m J ^ ^ ^ n L m 2m J 



p'^Rsip) 



-dp, 



^o(p) dp 



L m 2m J 

where Rs{p) is the radial charmonium wave function 
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As a result of the described transformations, we obtain the following general structure of 
pair Tjc production amplitude (|5]): 

M[gg r]^r]^] = ^(£1^2) + A2ieiP){e2P) + A3ieiQ){E2P) + 

A4{e,P){e2Q)+A5{e,Q){e2Q), ^ ' 

where Ai are the functions of variables s and t. Due to the bulkiness of corresponding 
expressions for coefficient functions Ai we do not present them here in exact form. 

In order to find the differential cross section for the gluonic subprocess we should calculate 
the squared modulus of the amplitude (fT3l) summed over polarizations of the initial gluons 
by means of the following relation: 

Then we obtain the general form of the gg — > rjcrjc cross section corresponding to the pro- 
duction amplitude (fT51) : 

^[99 ^ VcVc] = 1024 (^'-^1 + Ml + (^2 -A3-A, + A) {st + (M^ - tf)]'^ . 

(15) 

Making the substitutions for the functions Ai, we find it useful to transform the result ( !T5|) 
as follows: 

^[99^VcVc]{s,t) = ^^|i?(0)rX^u;,FW(.,t). (16) 

The auxiliary functions F^^^ entering the cross section ( |T6i) are written explicitly in Appendix. 
Note that the function F^^^ describes non-relativistic result which coincides in the limit 



M„ = 2m with the corresponding function obtained in Ref. [ll[ in the approach of NRQCD. 



Relativistic corrections in ( !T6|) are determined by a number of relativistic parameters Wj: 

h h 2 

Wo = 1, Wi = — , 0J2 = -r, UJ3 = UJi, 
-'0 -'0 



1,2 



, m + e{p) 2 J T fm + e{p) f m - e{p)\ . 2j 

= ' ^^^'P*" = J -MpT I^TTTm j (17) 





00 







In the non-relativistic limit, the parameter R{0) coincides with the definition of radial wave 
function at the origin, so it can be considered in some way as its relativistic generalization. 

All parameters, which contain the meson wave functions and describe the transition of the 
pairs (cc) to the bound state, are calculated in the framework of relativistic quark model [13, 



18| . This model is based on the Schrodinger equation with the Breit Hamiltonian in QCD 



and the nonperturbative confinement terms. Using the program of numerical solution of 



the Schrodinger equation [15|, |23|, we obtain relativistic wave functions and bound state 

fth 



energies of 5- wave charmonia. Numerical values of charmonium masses Mjy , = 3.072 GeV 
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TABLE I: The comparison of relativistic and nonrelativistic cross sections of a pair 5-wave char- 
monium production in pp collisions obtained for different sets of partonic distribution functions. 



Energy ^/S 


Meson pair, cross section type 


(T(total), nb 


a{2 < yp^ 


3 < 4.5), nb 






CTEQ5L 


CTEQ6L1 


CTEQ5L 


CTEQ6L1 


VS = 7 TeV 


J/iIjJ/iIj, relativistic 


9.6 


7.4 


1.6 


1.2 




J/il)J/'ip, nonrelativistic 


23.0 


17.7 


3.8 


2.9 




77c?7c5 relativistic 


23.7 


19.9 


1.3 


1.0 




r]cf]c-, nonrelativistic 


56.3 


48.1 


1.5 


1.2 


VS = U TeV 


J/ijjJ/ilj^ relativistic 


17.1 


13.2 


3.0 


2.1 




J/'il)J/'ip, nonrelativistic 


41.0 


31.6 


7.1 


5.1 




rjcTjc-, relativistic 


47.8 


39.3 


2.4 


1.7 




r/c^yc! nonrelativistic 


116.5 


94.7 


2.8 


2.0 



and M*|* = 2.988 GeV obtained in our numerical calculation lie close to the experimental 
results M't'jP = 3.097 GeV and M^'^p = 2.981 GeV M- The additional details on our 

J/lp Vc 



relativistic quark model can be found in Refs. pJl, Il8| . 

The numerical results of our calculation of the pair 5- wave charmonium production cross 
sections in the case of non-relativistic approximation as well as with the account of relativistic 
corrections of order v"^ are presented in Table [B Along with total cross section values, we 
have also included there the cross section predictions corresponding to the rapidity interval 
2 < yp,Q < 4.5 of the LHCb experiment [l| calculated with two different sets of linear PDFs: 
CTEQ5L jlll and CTEQ6L1 |26[. As shown in Table I, the cross section a[pp — )■ 2ric + X] 
at VS = 7(14) TeV is equal to 1.3 (2.4) nb for CTEQ5L and 1.0 (1.7) nb for CTEQ6L1. 
The most important production rates lie in the region of small Pt (see Fig. ([3])), where the 
color singlet contribution is dominant. Performing the numerical integration of differential 
cross section ( !T6|) . we use the LO expression for the running coupling constant as{fi) with 
the initial value as{Mz) = 0.118 and the renormalization scale fi = mx = 
where M is the meson mass. In our numerical calculations of the cross sections we set 
M^^ = 2.980 GeV and m = 1.55 GeV. Therefore, we take into account non-zero bound state 
energy of rjc charmonium state leading to the bound state corrections to the production cross 
section (fT6!) . Numerical results in Table I are determined by a number of parameters and 
functions: the c-quark mass, the factorization scale /i, parameters of the quark interaction 
operator, the bound state wave functions, the parton distribution functions and the strong 
coupling constant. Some of them (the c-quark mass, the quark- antiquark potential) are fixed 
in the relativistic quark model in the mass spectrum calculation. The factorization scale /i 
is taken in a commonly used form 1^, ll[ . Other quantities lead to basic uncertainties of 
our numerical results. 

It is evident from Table [H that relativistic corrections of order f ^ decrease the cross 
section values more than two times in both cases connected with the pair production of J/ijj 
or Tjc mesons. The only exception is the case of a pair rjc production cross section in the 
rapidity region 2 < yp^Q < 4.5, where the relativistic effects decrease the cross section only 
by approximately 20 percents. The change of PDF from CTEQ5L to CTEQ6L1 brings the 
additional 20-i-30 percent decreasing to the value of cross section. Along with the possibility 
of different PDF choices, there also exists the uncertainty dealt with the determination of 
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Pt, GeV 

FIG. 3: The cross section da/dPT at ^fS = 7 TeV for gluon distributfon functions CTEQ5L 
(dashed curve) and CTEQ6L1 (soHd curve). 



every particular partonic distribution function. The sets CTEQ5L and CTEQ6L1 contain 
no means to estimate the uncertainties of such sort, however the set CTEQ6M [26] has all 
necessary functionality. Using 40 uncertainty eigenvectors from CTEQ6M we can roughly 
estimate the error of every cross section value in Table [T] dealt with the PDF uncertainty 
as 15 %. The only known calculation of the ?7c-pair production in proton-proton collision 
was performed in Ref. ^if. Their Table II contains the obtained numerical results for 
different PDFs CTEQ5L and CTEQ6L1 with > 3 GeV, which are of order of 4 nb. 
They used almost the same values of the c-quark mass, the factorization scale as in our 
calculation but a different numerical value for the parameter i?(0). Our value of the radial 
wave function at the origin in nonrelativistic limit is equal to -R(O) ~ 0.8 GeV^^^, whereas 
in Ref. [ll[ the authors took the long distance matrix element {Oi)s = 1.4 GeV^, which 
gave R{0) = 1.7 GeV3/2. In the region of large transverse momentum Pt {Pt > 3 GeV) 
the cross section falls considerably, so that the basic contribution to our result in Table I 
is determined by the region of small momenta Pt- Therefore, our nonrelativistic results 1.5 
nb and 1.2 nb differ significantly from the values of cross sections obtained in nonrelativistic 
SPS approximation in |lli|. This difference is related with a choice of the parameter -R(O) 



m 



111 ] which exceeds our value more than two times. 
Another possible source of uncertainties is connected with the determination of relativistic 
wave function in the momentum region m. The obtained charmonium wave function is 
strongly decreasing in this region. Its numerical value at p = m is more than 50 times smaller 
the maximum value. Nevertheless, relativistic factors p"^ and p^, entering in the integrals 
/i^2 change this relation and increase the inaccuracy in the wave function determination at 
p ^ m. In spite of the fact that momentum integrals appear to be fully convergent, our 
relativistic model cannot provide a reliable calculation of the wave functions in the region of 
relativistic momenta p ^ m. Our definitions f[T7|) of the parameters Ji 2 contain the cutoff at 
relativistic momentum of order m. Using indirect arguments related with the mass spectrum 
calculation accuracy we estimate in 10% the uncertainty of the wave function determination. 
Larger value of the error would lead to the essential discrepancy between the experiment and 
theory in the calculation of the charmonium mass spectrum. Then the corresponding error 
in the cross section (flGjl is not exceeding 20%. We do not consider a part of theoretical error 
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related with radiative corrections of order because these corrections are omitted in our 
analysis. We also assume that relativistic corrections of order 0(f ^) to the cross section (fT6|l 
coming from the production amphtude should not exceed 30% of the obtained relativistic 
result. So, our total theoretical error is not exceeding 39%. To obtain this estimate we add 
the above mentioned uncertainties in quadrature. 
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Appendix A: The coefficients F^*) entering the differential cross section (I16|) 

In this appendix, we present analytical results for the parton differential cross section (fT6|l . 
Firstly, we introduce the following auxiliary functions of the Mandelstam variables s, t, u, 
and K = m/M: 



32 

?,M^sHu[2s + M2(l - 4^2)) (2t - M2(l + AK^)f [2u - M2(l + AK^)f ' 

-512 

kb = — — — : — Z-7T-, —. Z-7T-, — : — > (Al) 



9M^sHu{2s + M2(l - Ak^)) {2t - M2(l + Ak^)) {2u - M2(l + Ak^)) 

256{s-t + uf{s + t-u)^ 
27M^sHu{2s + M2 (1 - 4/€2)) ^ (2m - M2 (1 + AK^)f {2t-M^{l + Ak^)) 



6 ' 



ai = 3s^{t + u) + lASshu - 2s^{t + M)(3t^ - %%tu + 3^^) - 16s^t(19t^ - Mtu + 19m^)m+ 

/(t - uf(t + n)(3t^ - AAQtu + 3^^) + AsH(t - 5u){t - uf{5t - u)u + 212sH{t - u^x 

{t + u)u + 72st{t - u)^u + {k- 1/2) {9s^{t + u) + /(21t2 + 838tu + 21n^) + 2s^x 

(t + u){3t^ + 790tM + 3m^) - 2s^(9t^ + 5At\ - llA2t\'^ + 5Atu^ + 9u^) - 5s^x 

{t + u)(3t^ + 3A0t\ - A9At^u^ + 340^^^ + 3u^) - s^{t - uf{3t'' + gSOt^u + 626^^^^ 

980tu^ + 3u^) - l^sHit - ufiYIt^ - 86tM + \lu^){t + u)u - Astit - ufx 

(55^2 + 3Uu + ^'bu^)u - 72t{t - uf{t + u)u) , 

(A2) 



Aa = -165^^(125^ + 14/(t + u)- s^{30t'^ - 56tu + 30u^) - 27s^{t - uf{t + u) + 
As{t - 2u){t - uf{2t -u) + 9(t - u)^{t + u)) - 8stu{K - 1/2) (l43s^ + 245s^(t + u)- 
12s\5t^ - 34tu + Sn^) - 2s\t + u){109t'^ - 158tM + 109^^) - s^{t - ufx 
{m^ + 86tu + 17u^) + s{t - uf{t + M)(25t^ + 5Atu + 25u^) - lA{t - u)\t + uf), 

(A3) 
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6i = 3s^^{t + u) + 2380s^Hu - 12s^^{t + u){t'^ - 383tu + u'^) - As^Hux 

{U7lf - 5862tu + U71u^) + 2s^°(t + u){9t'^ - 807At^u + ISOSOi^i^^ - 807Atu^ + 9w^)+ 

4:sH{t - M)2(65t2 - WSUtu + 65u'^)u - 4s^(i - uf{t + u){3t'^ - U7m^u + 11992^^-^2- 

4479^1*^ + 3w^) + 20sH{t - u^iAm'^ + 1518tu + 4:13u'^)u + - u)'^{t + u)x 

{3t^ - 5312t^u + 22906*^1^2 - 5312tu^ + 3u^) - 32sH{t - uf{18Af + 251tu + 184^^)1^- 

8sH{t - uf{t + m) (259*2 - 22tu + 259u^)u + 8sH{t - uf{73t^ + 102tu + 73u'^)u+ 

1008s2*(t - + u)u + 288st{t - - {n - 1/2) (9s^^(t + u) - s^^x 

(39*2 ^ 33074i^i + 39^t2) _ 4^12^^ ^ m)(21*2 + 21470*ti + 21it2) + 4s^^(27*^ + 2386*=^w- 

85150*2x^2 + 2386*«=* + 27^x^) + 2s^°(* + ii)(99*^ + 99076*=^w - 270430*2^2+ 

99076*w=* + 99ii^) - 2s^(45*^ - 72738*^ - 101909*V + 414228*V - 101909*2-^^- 

72738*ii^ + Af>u^) - As^{t -uf{t + m)(45*^ + 18070*=^^ - 122638*2x^2 + I8070*w=^ + 45w^)+ 

4s^(* - 'u)2(3*^ - 34040*^^ + 26081*V + 83112*V + 26081*2^^ - 34040*^^ + 3m^) + 

s^(* - uY{t + 'u)(57*^ - 70124*^^ - 72506*2^2 _ 70124*^^ + 57m^) + s^(* - ufx 

{9f - 2350/' ?i - 63273*^^2 + 2204:t^u^ - 63273*2^^ - 2350*m^ + 9u^) + 24s^*(* - u^x 

(* + m)(1163*2 - 1906*M + 1163m2)m + 32s=^*(* - uf{562t'^ - 577f\ - 18*2m2_ 

577*^3 + 562m^)m + 40s2*(* - uf{t + m)(133*2 - 34*m + 133u'^)u + 144s*(* - u^^x 

(19*2 + mu + 19u'^)u + 864*(* - uY^it + u)u) , 

(A4) 

Ab = -4s2*m(1061s^^ + 1859s^°(* + u) - 6s^ (433*2 - 1502*n + 433^2) - 2s\t + u)x 
(2827*2 - 6134*M + 2827^2) + 24/(* - uf{72t^ - mtu + 72^2) + 40s^(* - uf{t + u)x 
(161*2 - 398*M + 161^2) + 2s^(* - m)^(313*2 + 5986*?;, + 313?72) - 2/(* - uf{t + u)x 
(1525*2 - 4586*M + 1525^2) - s^{t - m)6(1 189*2 + 3502*?;, + 1189^2) + s^{t- uf{t + u)x 
(261*2 - 1514*n + 261^2) + 12s(* - m)*(31*2 + 14*m + 31?72) + 144(* - uf^(t + u))- 
Astu{K- 1/2) (144215^2 ^342485ii(t + y) _ si0(8511*2 - 129586*?;, + 851 1^2)- 
8s^(* + m)(9129*2 - 23527*M + 9129^2) - 2s®(16065*^ + 59372*3^ - 167130*2^2+ 
59372*^3 + 16065m^) + As\t - uf{t + m)(10827*2 - 45008*?;, + 10827?i2) + 2s^(* - ufx 
(18891*^ + 2840*3^ - 77062*2^2 + 2840*^3 + 18891m^) + 52s^(* - uf{t + u)x 
(97*2 ^ ^2{)tu + 97^2) - s4(* - m)4(6659*^ - 20876*3m - 3822*2^2 - 20876*?/,''^ + 6659?/)- 
4s3(* - vf{t + ?7,)(1869*2 - 2300*?x + 1869?/) - /(t - uf{3l9lt* + 1344*-\ - 9454*2?x2+ 

1344*^3 + 3191m'^) - 4s(* - uf{t + -u) (203*2 + 334*m + 203^2) - 432(* - uf^{t + uf) , 

(A5) 
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ci = 3si^(t + u) + 127924:S^Hu - 4:S^^{t + u){3f - 99607tM + Sm^) - 16s"tx 
(6327t2 - 119146tM + 6327u^)u + 2s^'^{t + M)(9t^ - 820440^^^ + 1900702tV- 
820440tM^ + 9u^) - 32s^t(42787t^ + 95828t^'u - 266062t2M2 ^ 95828^^^ + 42787m^)m- 
4s^(t - + M)(3t^ - 476492^3^ + 1778482tV - 476492^^^ + 3^4^, ^ g^r^j^^ _ ^^,2^ 
(395023t^ + 173924t\ - 1092582^21*2 + 173924tM2 + 395023m^)m + s%t - u^x 
(t + M)(3t^ - 45892^3^ + 6221826tV - 45892^^3 + 3^4^^ _ 4^5^^^^ _ ^^,4^ 
(587629t^ - 129076t^M - 999794^^^2 _ i29076tM^ + 587629m^)m - Ash{t - ufit + u)x 
(258537^2 + 588302tM + 258537m2)m + 8sh{t - M)^(56911t^ - 55796t\ - 45238^2^2- 
55796tM^ + 56911m^)m + 1440s2t(t - + M)(287t2 + 166tM + 287^2) + 

288st(t - M)^%(255t2 + 538tM + 255^2) - 864t(t - u^^it + u)u, 

(A6) 

Ac = -8tu(26215s^^ + 72029s^^{t + u) - s"(46586t2 - 347812tM + 46586^2)- 
4s^°(t + M)(78089t2 - 168598tu + 78089^2) - 4s^(37679t^ + 168700t\ - 401590t2u2+ 
168700tu3 + 37679u^) + 2s^{t - uf{t + u){223529t'^ - 615866tu + 223529^2) + 
2s\t - u)2(227477t^ + 195492t\ - 823282t2u2 + 195492tu^ + 227477u^)- 
4/(t - u)\t + u)(50509t2 - 265278tu + 50509^2) - s\t - M)^(402183t^ - 13884t^u- 
817942t2u2 - 13884tu^ + 402183u^) - 7s\t - u)%t + u) (6701^2 + 56246tM + 6701^2) + 
24s^(t - M)^(4913t^ - 2740t^M - 6138^2^2 - 2740tu^ + 4913u^) + 48s2(t - uf{t + u)x 
(817^2 + 812tu + 8171/2) - 324s(t - n)^°(llt2 + 2tu + llu'^) - 972(t - u)^^{t + u)), 

(A7) 

02 = ai + Aa, 1)2 = hi + Ah, C2 = Ci + Ac. (A8) 

Then, the coefficients F^^^ entering the differential cross section flTB]) have the following 
form: 

FW = A:2(a? + a2), 
F'-^^ = kahiaihi + a2&2) + 8A;2(a2 + aj), 

F(2) = -4A;2(a? + a2), (A9) 

= Qkah{aihi + 02^2) + kakciaici + a2C2) + 24A;2(a2 + a^) + ^kl{hl + hi). 

Note that we expand functions ai^2 and 61,2 in the mass difference (2m — M) up to the 
term linear in {k — 1/2) and set the value k = 1/2 in ci,2- Such simplifications allow us 
to significantly reduce the length of analytical expressions f lA2p - flA7p . while the numerical 
results of the cross sections change on 1 ^5 percents. In Table [T] we present numerical results 
corresponding to exact functions Oj, hi, and Cj. 
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